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Abstract 

We investigate a model of interacting clusters which compete for growth. For 
a finite assembly of coupled clusters, the largest one always wins, so that all but 
this one die out in a finite time. This scenario of 'survival of the biggest' still holds 
in the mean-field limit, where the model exhibits glassy dynamics, with two well 
separated time scales, corresponding to individual and collective behaviour. The 
survival probability of a cluster eventually falls off according to the universal law 
(lnt)" 1 / 2 . Beyond mean field, the dynamics exhibits both aging and metastability, 
with a finite fraction of the clusters surviving forever and forming a non-trivial 
spatial pattern. 
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1 Introduction 



Non-equilibrium dynamics can lead to counter-intuitive situations. One is used, for 
example, to the premise of equilibration: in a system with unequally distributed masses, 
the effect of most 'physical' interactions would be to bring the system to an equilibrium 
state where masses are distributed equally. This equilibration principle is known to fail 
in some physical instances, mostly in the presence of long-range forces, the prototypical 
example being gravitational forces. In fact, as is well known, the effect of gravitation is 
to amplify forever the contrasts in mass distribution throughout the Universe £Q. 

In this work, we interest ourselves in an extreme case of disequilibration. The model 
investigated below deals with immobile interacting clusters which compete for growth. 
Although it arose in an astrophysical context, that of mass accretion by black holes 
coupled by the radiation field in a brane world [21 E] , its emergent features are relevant 
to a far wider range of problems. The present model is strongly out-of-equilibrium 
and it obeys no mass conservation law. A variety of transient behaviour is therefore 
possible: for example, two interacting clusters can both decay, or both grow before one 
of them dies out. At late stages, the model follows the survival of the biggest scenario, 
an example of Darwinism in a physical system. In the mean-field geometry, the largest 
cluster generically wins out over all the rest; in finite dimensions, one has the possibility 
that infinitely many clusters survive and grow forever, provided each of them is isolated, 
in a sense that will be clearer later on. This is actually quite meaningful in the original 
astrophysical context, since it could, with appropriate modifications, be adopted to model 
the scenario of primordial black holes evolving to a size such that they survive in the 
present era. 

The present model, defined in Section [21 deals with an assembly of pointlike, im- 
mobile clusters, which are entirely characterised by their masses. Cluster masses evolve 
according to coupled deterministic, non-linear first-order equations. We address a med- 
ley of situations ranging from finite assemblies of clusters, to the thermodynamic limit, 
examined both in the mean-field geometry and on a lattice with nearest-neighbour in- 
teractions. Section [3] describes the dynamical behaviour of a single isolated cluster: a 
large enough cluster, whose initial mass exceeds some threshold, grows forever, whereas 
a smaller one evaporates and disappears in a finite time. Section 0] concerns our find- 
ings on two interacting clusters, and more generally finitely many coupled clusters: the 
generic scenario is then the survival of the biggest, so that the largest cluster wins out 
over all the rest. The mean-field regime of a large collection of weakly coupled clusters 
is investigated in Section |3| The system exhibits aging and glassy dynamics, involving 
two well-separated time scales. The cluster survival probability decays according to the 
universal law (hit) -1 / 2 . In Section El we examine the model with nearest-neighbour in- 
teractions in finite dimension: the dynamics now exhibits both aging and metastability. 
The finite fraction of survivors, i.e., clusters which survive and grow forever, builds a non- 
trivial spatial pattern. In the Discussion (Section |7|), we put our results in perspective 
with other growth models. 
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2 The model 



The model investigated in this work is a direct generalisation of that derived in [21 E|- 
Consider n pointlike, immobile clusters, which are entirely characterised by their time- 
dependent masses mj(i) for i = l,...,n. The cluster masses evolve according to the 
following coupled deterministic, first-order equations: 

dm; I a 1 ^ dmA 1 

These dynamical equations were originally written to model the kinetics of black hole 
growth in a radiation fluid [21 El- In that context, they only hold after some microscopic 
initial time to- The positive (gain) term in the right-hand side of (j2.1|) represents mass 
accretion by the black hole from the surrounding fluid. The accretion rate in the large 
parenthesis is the sum of the free rate for an isolated black hole, proportional to the pa- 
rameter a > 1/2, and of the rate induced by all the other black holes via the surrounding 
fluid. The coupling between black holes i and j is proportional to the inverse square 
distance <i^(to) between them at the initial time to- The negative (loss) term in the 
right-hand side of (J2.1)) represents evaporation due to Hawking radiation. 

In this paper, the dynamical equations (j2.1|) . or equivalently ()2.4|) . are now seen as 
representing competitive cluster growth. The functional form of the original equations 
derived in |3] is kept unchanged for definiteness. The symmetric matrix of couplings g^ 
reflects the underlying geometry. Succeeding sections will deal with a raft of scenarios of 
ever-increasing complexity, ranging from the dynamics of two coupled clusters to that of 
infinitely many. 

It turns out to be convenient to switch from physical time t to reduced (logarithmic) 
time 

s = lnf, (2.2) 
to 

so that the initial time to is mapped onto the origin s = 0. Furthermore, we introduce 
for convenience the reduced masses and square masses: 

m i 2 _ m l ( o\ 

' ~~ £l/2' Di — X i — ^ ■ 

The dynamical equations (j2.1|) then become the following autonomous equations 
dx; , I 2a — 1 ^ fXi ,\\ 1 



for the reduced masses Xi(s), which exhibit no explicit dependence on the reduced 
time s [3]. Throughout the following, accents will denote differentiation with respect 
to the reduced time s. It will also be assumed that the couplings #y are small enough, 
so that 

det (Sij + g i:j .r, ) ; , j | „ > 0. (2.5) 
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When this inequality holds, the time derivatives x\ can be solved explicitly from the im- 
plicit dynamical equations (J2.4j) . so that the dynamics is regular. Whenever the regularity 
condition (|2.5|) holds at s = 0, it turns out to be preserved by the dynamics. 



3 One isolated cluster 

The simplest situation is that of a single isolated cluster of mass m(t). The dynamical 
equation (|2.1|) reads 

dm am 1 

~TT = ~4- ' (3.1) 

at t m 

The dynamical equation (|2.4j) simplify to the following ones for the reduced mass x(s) 
and square mass y(s): 

. 2a - 1 1 . . 

x'=—— x--, 3.2 
2 x 

y' = (2a - l)y - 2. (3.3) 
Equation ()3.3|) is the easier to solve. It yields at once 

y (s) = y, + (y -yM 2a ' l) \ (3-4) 

where m = m(t ) and y = m^/t are the initial values of m(t) and y(t), respectively, 
whereas 

is the unstable fixed point of ()3.3|) . 

Returning to physical variables, (|3.4|) reads 

/ t \ 2a 

m{tf = y4 + {m 2 - yj ) i^-J . (3.6) 

This implies that for all a > 1/2, we have two kinds of behaviour: 

• Large clusters, whose initial mass is such that yo > y±, i.e., mo is larger than the 
mass threshold 

^ = W 1/2 =fA) 1/2 . (3-7) 



,2a - 1, 

are survivors: they survive and keep on growing forever. Equation (j3.6J) reads 
alternatively 

t f t \ 2a 

m(t) 2 = ml- + (ml-ml)(-) . (3.8) 
to \to/ 

The second term is the leading one at late times, for all a > 1/2. 
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• Small clusters, whose initial mass is below the threshold: y < y*, i.e., m < to*, 
evaporate and die out in a finite reduced time, 

s(y ) = -^—\n^—, (3.9) 
2a - 1 y* - y 

which diverges logarithmically as the mass threshold m± is approached from below. 
The corresponding physical time, 

/ \ V(2«-D 

tW = U— , (3.10) 

\y* -yoj 

diverges as a power law. 

Even in this simple case of independent clusters, we get an indication of a Darwinian 
scenario (where size, i.e., mass replaces fitness): the biggest clusters survive, while the 
smaller ones die out. 

Consider now a very large assembly of isolated, i.e., non- interacting, clusters, charac- 
terised by the (continuous) probability distribution function P(yo) of their initial square 
masses yo. One of the quantities of most interest is the survival probability S(s), defined 
as the fraction of the clusters which have survived up to reduced time s. Surviving clusters 
are exactly those whose initial square mass obeys y > Y(s), where the time-dependent 
threshold Y(s) is the inverse of s(yo) introduced in ()3.9|) : 

Y(s) = (l - e-* 2 *" 1 )') y*. (3.11) 

The survival probability at time s reads therefore 

POO 

S(s)= P(y )dy . (3.12) 

JY(s) 

The limit survival probability 1 is defined as the fraction of survivors, i.e., clusters 
which survive and grow forever. These are the clusters whose initial mass is above the 
threshold y+ introduced in (|3.7|) . We thus obtain 

POD 

S {1) = lim S(s) = / P(y )dy . (3.13) 

For simplicity, we shall often consider in the following an exponential distribution of 
initial square masses: 

P(l/o)=/ie-w». (3.14) 

The result (I3.13J) then reads 

S' (1) =e-«'*. (3.15) 



1 The subscript (1) recalls that this result holds for isolated clusters. 



5 



4 Two interacting clusters 



Until now, we have considered only independent clusters: these survive or not, depend- 
ing on their initial masses. We now turn to the more interesting situation of interacting 
clusters. The form of the interactions was derived in earlier work [3| and is such that 
clusters could 'feed on' each other: thus, smaller clusters disappear faster as if they were 
swallowed by the larger ones. In this section, we explore the details of the simplest pos- 
sible case, that of two interacting clusters with masses m\(t) and rri2(t), and interaction 
strength g = g 12 = g 21 . 

We look successively at the special case of equal masses (Section 4.1) and at the 
generic case of unequal masses (Section 4.2). Only in the first case are the two masses 
able to survive forever, growing more slowly than if they had been alone. In the second 
case, the bigger cluster swallows the smaller one, generically. 

The dynamical equations ()2.4|) read 



, ( 2a -I (x 2 ,\\ 1 

x i = — n 9 hr + x 2 )) x i 



V 2 3 V2 ' VJ L Xl 
. (2a -1 (x x ,\\ 1 

Solving these equations for the time derivatives, we obtain 

(2a — l)x\x 2 — 2x 2 + 2g(l — ax 2 2 )x\ + g 2 x\x\ 

2x 1 x 2 (l - g 2 x 1 x 2 ) 
{2a — 1)0:1X2 — 2x\ + 2g(l — ax\)x\ + g 2 x\x\ 



(4.1) 



(4.2) 



,X 9 



2x x x 2 {\ -g 2 x 1 x 2 ) 



In the generic situation of two unequal masses, the above dynamical equations already 
illustrate the full complexity of the problem. The regularity condition ()2.5|) . which reads 

1 - g 2 x x x 2 > (4.3) 

in the case of two clusters, is needed for the denominators not to vanish. 



4.1 Equal masses 

The main results of this section can be explained physically in the following way: since 
the interactions cause each mass to 'feed on' the other, overly strong interactions will 
lead to a strongly depletive effect on both, as a result of which neither survives. On the 
other hand, a weakly interacting pair of equal mass clusters can, provided their masses 
are above a threshold, survive in gentle symbiosis; both depletion and accretion keep 
occurring at comparable rates, and the pair survive forever. 

Consider now two clusters whose masses are equal at the initial time t . This symme- 
try is clearly preserved by the dynamics. Let x(s) be the common value of their reduced 
mass. Equations ()4.2j) simplify to 

(2a - l)x 2 -2-gx* 

2x(l + gx) y ' 1 
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The fixed points of the above dynamical equation, given by 



(2a - l)x 2 - 2 - gx 3 = 0, (4.5) 

dictate the qualitative features of the dynamics. There is a critical value of the coupling, 

/2(2a- 1) 3 \ V2 , x 

9c = 1 n „ ' , (4.6) 
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which separates two kinds of behaviour: 



• For a large enough coupling (g > g c ), ()4.5|) has no real positive root, leading to the 
absence of a physical fixed point. The reduced mass x(s) of both clusters decreases 
monotonically until they simultaneously die out in a finite reduced time. 

• For a small enough coupling (g < g c ), ()4.5|) has three real roots, two of which are 
positive and correspond to physical fixed points: 

y l J 2 < X(i) (unstable) < (3y*) 1/2 < x@) (stable). (4.7) 

Small clusters, such that x < xm, are attracted by x = 0, so that both disappear 
in a finite time. Large clusters, such that x > X(i), are attracted by x^y. those 
pairs of clusters are survivors, and their common mass grows as 

m(t) « x (2) t 1/2 . (4.8) 

Note that this growth rate is slower than that of an isolated cluster [see (13. 6|) ]. 

In other words, small masses, no matter what the coupling strength, die out in a finite 
time. For large reduced masses, the role of the coupling strength matters. For g < g c , 
those larger than xn\ survive forever, and grow more slowly than if they had been isolated; 
for g > g c , all die eventually. 

These results match those of Section|3]in the limit of a vanishingly small coupling. The 
unstable fixed point has a finite limit xm — > yV 2 , whereas the stable one diverges as xm ~ 
(2a — l)/g. We recall from Section 01 that y± is the threshold above which independent 
clusters survive forever. The unstable fixed point xm, which likewise separates dying 
from surviving clusters, thus reduces to y* in the g — > limit, as it should. The main 
new element with respect to the non-interacting limit is the existence of the stable fixed 
point X(2)- However, its effect is to attract all masses above Xm to itself, so that here, 
too, the effective behaviour is unchanged with respect to the noninteracting case; masses 
above xm survive forever. 

Finally, it is worth coming back to the regularity condition (|4.3J) . In the present 
situation of two equal masses, this condition reads 

X < 34im) -^lim — • (4-9) 

9 
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The limiting value x\\ m plays no special role in the dynamics of two clusters with equal 
masses. It will however play an important part in the transient dynamical behaviour of 
two unequal masses, to be studied in Section 4.2. It is therefore worth investigating the 
fate of xi im . Equation (|4.4|) implies 



x 



X X\\i 



a-1- g 2 
2g 



This expression singles out the following value of the coupling strength: 

l) 1 / 2 (a > 1). 



9\h 



a 



(4.10) 



(4.11) 



For g < gu m , the right-hand side of ()4.1(Jj) is positive, so that x\ im flows toward larger 
values of x. Conversely, for g > g\ im , the right-hand side of ()4.1Uj) is negative, so that x\\ m 
flows toward smaller values of x. We present in Figure ^ for further reference, the 
phase diagram of the two-cluster problem in the a-g plane. Four different phases can be 
defined (see caption), according to the number of real positive fixed points, and to the 
number of those obeying the condition (|4.9J1 . The phase boundaries are determined by 
equations (JUj) and (l4~TTl) . 
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Figure 1: Phase diagram of the two-cluster problem in the a-g plane. Full line: critical 
coupling g = g c (a) of (|4.6j) . Dashed line: g = gn m (a) of (|4.11|) . Phase I: no real positive fixed 
point. Phase II: two fixed points, both obeying ()4.9|) : Xm < < x y m . Phase III: two fixed 
points, only xn\ obeys (|4.9|) : xn\ < xn m < xm. Phase IV: two fixed points, none obeying (|4.9jl : 
x lim < < X(2)- Full symbol: quadruple point (a = 5/4, g = 1/2): xn\ = xr 2 ) = ^lim = 2. 
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4.2 Unequal masses 



In the general case where the clusters have unequal masses, the role of the interactions 
is inherently disequilibrating: mass differences, however small initially, get rapidly am- 
plified, leading to the generic scenario of the survival of the biggest. We shall investigate 
successively three stages in the dynamics of two clusters with slightly unequal masses. 



Early stage: linear stability analysis 

Consider first the early stage of the dynamics for two clusters with a small mass difference. 
Setting 

xi(s) = x(s) + e(s), x 2 (s) = x(s) — e(s), (4-12) 

to first order in the difference e(s), the mean reduced mass x(s) obeys (|4.4j) . while the 
difference itself obeys the linear equation 

e'(s) =L(x(s))e(s), (4.13) 

where the instantaneous Lyapunov exponent h(x) reads 

L( X ) = U(,) + e{2 r t 9X+ 2 a f\ U W = ^ + i (4.14) 
x(l — g l x A ) 2 x z 

The full Lyapunov exponent L(x) of two interacting clusters is therefore larger than the 
Lyapunov exponent L (x) in the absence of coupling, which is in turn larger than the 
constant (2a — l)/2. As underlined above, interactions thus enhance disequilibration. In 
any case, irrespective of the mean initial mass and of the coupling, any small initial mass 
difference diverges exponentially in the early stages of the dynamics. In particular, the 
fixed point x(%) of Section 14.11 which is stable against a symmetric perturbation of the 
form Sx\ = 8x2, is always linearly unstable against an asymmetric perturbation of the 
form Sxi = —5x2 = £■ 



Intermediate stage: transient behaviour in the various phases 

Later stages of the dynamics cannot be described in closed form, because of the non- 
linearity of (|4.2j) . The detailed transient time dependence of both masses depends on the 
location of the parameters a and g in the phase diagram of Figure ^ especially when the 
initial mass difference is small. 

These features are illustrated in Figure showing the shape of typical trajectories 
in the gx\-gx 2 plane. The dashed line shows the limit of the regularity condition ()4.3|) . 
so that allowed pairs of reduced masses are below that line. Each full line shows a 
trajectory starting with a small mass difference e = ±10~ 2 . The mass difference then 
grows monotonically, until the trajectory hits either of the co-ordinate axes in some finite 
time si, when the lighter mass disappears. Each panel corresponds to a typical choice 
of a and g in each of the four phases. Phases I and IV are very similar: the lighter 
mass always decreases monotonically, whereas the larger one decreases in a first stage. 
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If the larger mass is large enough, it may then start increasing before the lighter dies 
out. Phases II and III are also similar: both masses may increase in a first stage if their 
difference is small enough. 




Figure 2: Typical trajectories in the gx\-gx2 plane, for one choice of the parameters a and g 
in each of the four phases shown in Figure ^ Full lines: a few trajectories with an initial mass 
difference e = ±10 -2 . Dashed line: limit of the regularity condition (|4.3|) . i.e., gx\ ■ gx2 = 1. 
Full symbols: fixed points. 



Late stage: survival of the biggest 

The values of the parameters a and g become asymptotically irrelevant in the late stages 
of the dynamical evolution. Indeed, as illustrated above, any trajectory eventually hits 
either co-ordinate axis. There is therefore one single generic scenario for two clusters 
with unequal initial masses, i.e., that of survival of the biggest: 

• The smaller one dies out in a finite time s±. 

• The larger one then evolves according to the results of Section EJ Depending on 
the value of its mass at reduced time s±, it may either also disappear in a finite 
time (for y(si) < j/*), or survive and grow forever (for y(sx) > y*). 
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The above results still hold qualitatively for any finite number n > 2 of fully inter- 
acting clusters (all couplings are non-zero). In the generic situation of unequal masses, 
the scenario of survival of the biggest applies: the n — 1 smaller clusters die out one after 
the other, while only the largest one may become a survivor. 

5 Mean-field limit 

Having explored the behaviour of finitely many interacting clusters, we now turn our 
attention to collective behaviour in the thermodynamical limit of an infinite assembly 
of interacting clusters. In the present section, we focus on the mean-field regime of 
long-range interactions. The main feature of mean-field dynamics is again that all the 
clusters eventually die out, except the largest one. The clusters which survive up to 
time s are those whose initial reduced square mass exceeds some time-dependent mass 
threshold Y(s), to be determined below. The general case is necessarily somewhat formal 
(Section 5.1). If, however, one considers the regime of weak interactions, the formalism 
simplifies considerably (Section 5.2). 

5.1 General formalism 

We consider the mean-field limit of a large assembly of clusters (n^> 1), assuming that all 
the couplings have the same value g. We perform the usual rescaling of the interaction 
strength in mean-field models: 

g = l. (5.1) 

n 

The problem simplifies drastically in the thermodynamic limit, defined as usual as 
the n — > oo limit at fixed g. In this limit, the coupling strength of any cluster to its 
whole environment, measured by g, remains of order unity, whereas the strength of the 
coupling between any two different clusters, measured by g, falls off as 1/n. 

In this thermodynamic limit, (|2.4jl implies the following dynamical equation 

y'(s) = 7(s)y(s) - 2 (5.2) 
for the reduced square mass y(s) of any of the clusters. We have introduced the notation 

7(s) = 2a - 1 - g (M(s) + 2M\s)) (5.3) 
for the effective growth rate of the square mass y(s), where 

M{s) = (x) s = (y 1 / 2 ) s = Jim i £ y^s) 1 ' 2 (5.4) 

is the mean reduced mass of the clusters at reduced time s. 

Despite its apparent simplicity, (|5.2|) is non-trivial, because of its self-consistency: 
its right-hand-side indeed involves j(s), and therefore M(s), and therefore the solution 
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of ()5.2j) itself. This self-consistent problem can be solved formally as follows. First, we 
have on differentiating (J5.4j) 

M'(s) = ^M(s)-N(s), (5.5) 

with 

N(s) = (y^ 2 ) s =\im±- Y, (5-6) 

i:j/j(s)>0 

where only non-zero values of y^, corresponding to clusters i which survive at time s, are 
involved in the sum. The effective rate 7(s) can therefore be solved from (J5.3|) : 

2a-l+g(2N(s)-M(s)) 

7(s) = TTm*) ' ( } 

On the other hand, a formal integration of ()5.2|) yields 

y(s) = G(s)(y -Y(s)), (5.8) 

with 

GW=«p(j£ •*.>*.), ^1 = 2/;^, (5.9) 

so that 

GW = j^. 7(.) = (5-10) 

These steps lead to the following picture of the mean-field dynamics. For a given 
initial value y = y(0) of the reduced square mass, the solution ()5.8|) holds only as 
long as y(s) is positive, or, equivalently, yo > Y(s). Hence the smaller clusters, with 
initial square masses yo < Y(s), have already disappeared at reduced time s, while larger 
ones, such that yo > Y(s), have their square masses shifted and dilated from y to y(s), 
according to ([5.8)1 . 

This has strong echoes of the case of many non-interacting clusters. Recall that there 
existed a mass threshold, also called Y(s) in ([3. lip , below which all particles had died at 
time s, and above which they survived. The quantity Y(s) in (|5.8p above generalises this 
threshold in the presence of a mean-field coupling. Now, clusters below this threshold die 
as before, while the mass of a cluster grows, as a result of interactions, from y to y(s). 
In the absence of coupling (g = 0), the present Y(s) reduces to that obtained in (|3.11jl . 

The above formalism allows us, provided the continuous probability distribution of 
initial square masses P{yo) is known, to express all the quantities of interest in terms of 
a single (so far unknown) dynamical quantity, the threshold Y(s). In terms of this, M(s) 
and N(s) read 

M(s) = G(s) 1 / 2 / (y - Y(s)f 2 P(y ) dy , (5.11) 

JY(s) 
roc 

N(s) = G(s)-^ 2 j Y Jyo~Y(s))-^ 2 P(yo)dyo. (5.12) 
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Equation (|5.7|) then provides a self-consistent equation for the unknown quantity Y(s), 
as it only involves Y(s) itself and its first and second derivatives. Of course, the resulting 
non-linear integro-differential equation cannot be solved in closed form in general. 

To recapitulate, the program for the mean-field solution of the dynamics of n ^> 1 
clusters is the following: Equation ()5.8|) gives the growth law of any cluster mass in 
terms of the time-dependent mass threshold Y(s). Only clusters above this threshold 
survive at reduced time s, as the rest have died. The calculation of Y(s) can be done 
self-consistently, at least in principle, from (|5.7j) . Quantities of interest, such as the 
mean mass of surviving clusters at time s, can then be calculated via (j5.11|) and similar 
expressions. 

The definition of the survival probability S(s) is unaltered by the presence of inter- 
actions, and is still given by equation ()3.12jl . i.e., 



where ((• • -)) s denotes the normalised mean over the clusters surviving at reduced time s. 
In general, average quantities computed over all clusters, dead and alive, must be renor- 
malised by S(s) in order to get an appropriately normalised average over surviving clus- 
ters at any time s. 

5.2 Weak-coupling regime 

The results of Section 5.1 hold for arbitrary values of g, and their formal nature admits 
of no further simplification. However, a much greater transparency is achieved in the 
regime where the rescaled coupling g is small. 

In this regime, the dynamics consists of two successive stages. In Stage I, the clus- 
ters behave as if they were isolated, i.e., their masses evolve according to the results of 
Section El This fast stage of the dynamics therefore corresponds to individual behaviour. 
The only surviving clusters after Stage I are those whose initial masses exceed the thresh- 
old (|3.7|) . The effect of interactions sets in at Stage II. This slow stage of the dynamics 
corresponds to collective behaviour. All but the largest cluster eventually die out during 
this stage. 

The weakly interacting mean-field regime of our model of interacting clusters therefore 
exhibits characteristic features of glassy systems 4j. These aging phenomena originate 
in the presence of two well-separated time scales of fast and slow dynamics, with a 
ratio of respective time scales growing as l/g 2 . Another striking feature of our model 
is the universality of the main asymptotic results in Stage II dynamics: the survival 
probability falls off generically as (hit) -1 / 2 , whereas the mean square mass of survivors 
grows as tint. As is well known such logarithmic behaviour is another telltale sign 
of glassy dynamics. We describe all of this below in more detail. 




(5.13) 



The mean reduced mass of the surviving clusters reads 



W. = ((y 1/2 ))s = 




(5.14) 
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Stage I: Fast individual dynamics 



In this first stage, interactions are essentially irrelevant, and the dynamics is fast. The 
mass of each cluster evolves according to Section El independently of all the others, as if 
it were isolated. The survival probability decays from 5(0) = 1 to the plateau value 
of ()3.13jl . whereas the time-dependent threshold Y(s) of (J3.11)l increases from Y(0) = 
to y* of dS3J). 



Stage II: Slow collective dynamics 

In this second stage, the interactions are responsible for a slow collective dynamics in the 
weak-coupling regime. 

The evolution throughout Stage II can be described as follows. Starting from the 
assumption (to be checked later on) that the dynamics is slow, we have 7(5) <C 1 and 
M'(s) <C M(s). Equation (J5.3)l therefore simplifies to 

. . 2a — 1 / 1 

M(s) ss . 5.15 

9 

Now, inserting f)5. lUj) and (J5.15)l in ([5.11)1 . we obtain the following closed differential 
equation for the unknown Y(s): 

where the function 

/CO 
(y -Y)^ 2 P(y )dy (5.17) 

is entirely determined by the initial mass distribution. 

The behaviour of the threshold Y(s) throughout Stage II is obtained by integrat- 
ing (J5.16)) . with an initial value equal to the plateau value y+ of (|3.5|) : 



y dy 2g 2 s 



R(yf (2a -If 

Equation 1)5.18)1 contains the key to the dynamical behaviour in Stage II. 
First, the characteristic time of the collective dynamics, 

(2a - l) 2 



(5.18) 



(5.19) 



becomes arbitrarily large in the weak-coupling regime (lj — > 0). The assumption of slow 
dynamics is thus fully justified. 

Then, the analysis of the long-time dynamics goes as follows. Note that the left-hand 
side of ()5.18)1 diverges as s — > 00. The only way this can occur is if R(Y) falls off to zero 
for long times, i.e., if the threshold Y(s) goes to the maximum possible initial mass ?/ max , 
i.e., more formally, the upper bound of the continuous distribution P(yo)- The survival 
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probability S(s) then also falls off to zero for long times. We can conclude that the whole 
population of clusters which survived Stage I eventually disappears during Stage II of 
the mean-field dynamics. For a finite mean-field system consisting of n coupled clusters, 
at most one of them will survive forever, as already mentioned at the end of Section 4.2. 

Having established the general pattern, we now specialise to specific distributions of 
initial masses, in order to obtain quantitative predictions. We first consider an exponen- 
tial distribution of initial masses. We will find that the results obtainable from it can be 
generalised to a raft of other distributions. 

With the exponential distribution (|3.14|) . equations (|5.13j) . (j5.18J) yield 

~ S(sY 6 + (2a -1)^ 
In the late times of Stage II, the survival probability therefore decays as 

S{s)^ 2 -^r±{Cs)-^\ (5.21) 

with 

C = n (5.22) 

for the chosen exponential distribution, irrespectively of a, fi, and g, provided the latter 
is small enough. 

Equations (PHD . (15~T5T) then lead to 

«*». « {Csf'\ (5.23) 

For a final presentation of the above results, we return to physical variables. In terms 
of these, the survival probability falls off as 

2<r — 1 / t \ - 1 / 2 
S(t) w — — [C In- , (5.24) 



9 V t / 

whereas the mean mass of the surviving clusters grows as 



t \ 1/2 

(m)) t w [Ct In- . (5.25) 



The universality inherent in the scaling results ()5.21|) - (j5.25|) is unusual, because it 
includes the prefactor C, which is itself independent of the details of the initial distribu- 
tion P(yo) of square masses. It can indeed be checked explicitly that C only depends on 
the tail exponent of this distribution in the vicinity of its upper bound y max : 

• In the bounded case (y max finite), assuming that the distribution has a power- law 
behaviour P(yo) ~ A(y max — yo) a l f° r Vo ~^ Vmax, with a tail exponent a > 0, we 
obtain 

C = na\^±^\ . (5.26) 



r ( a + l) 
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In the unbounded case (y max = oo), assuming that the distribution has a power- law- 
behaviour P(yo) ~ B y^ 6-1 for y — > oo, with a tail exponent b > 1/2, we obtain 

C = -M^7^l • (5-27) 

It is worth noticing that both results (|5.26|) and (|5.27j) smoothly converge to the 
particular value C = n of (|5.22jl . as the tail exponents a and b get large: 

C = »(l-^ + ...)=»(l + » + ...). (5.28) 



6 Finite-dimensional lattices 

The mean-field regime described above can be thought of as an infinite-dimensional limit 
of our model. Now, in order to include the effects of fluctuations, we consider a finite- 
dimensional lattice model. In addition to the emergence of two well-separated time 
scales, the model now displays metastability: the system gets finally trapped forever in a 
non-trivial attractor, where every surviving cluster is isolated. 

More specifically, clusters sit at the vertices n of a regular lattice. Every pair of nearest 
neighbours interacts with a uniform coupling strength g. Numerical simulations have 
been performed on hypercubic lattices: the chain (D = 1), the square lattice (D = 2), 
and the cubic lattice (D = 3). The coordination number of these lattices is z = 2D. 
Throughout this section, initial masses are given by the exponential distribution ([3. 14)1 . 
Unless otherwise stated, we set a = 1, g = 10~ 4 , and \i such that (|3.15|) yields = 0.9. 

The main focus will be the limit of weak coupling (g < 1). We therefore simplify the 
dynamical equations (|2.4|) . by keeping terms up to first order in g. The resulting explicit 
equations 

+ 9^2(— aim)]^--, (6.1) 

where m runs over the z nearest neighbours of site n, are solved numerically by means 
of a standard first-order scheme. 




6.1 Two-step dynamics 

In the weak-coupling regime, the dynamics generated by ()6.1|) again consists of two 
successive well-separated stages. As before, fast individual dynamics are exhibited in 
Stage I, while Stage II is the arena for slow collective dynamics. The effects of going 
beyond mean field are only palpable in the latter stage, since interactions are irrelevant 
in Stage I. 
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Stage I: Fast individual dynamics 

There is little new here with respect to the mean-field regime. The mass of each cluster 
evolves as if it were isolated, as before. The survival probability S(s) decays rather fast 
from S(0) = 1 to its plateau value of ()3.13|) . 



Stage II: Slow collective dynamics 

The slow collective dynamics throughout Stage II is now very different from the mean- 
field regime. The survival probability S(s) indeed decays from its plateau value Sm to a 
non-trivial limiting value S'(oo), because of metastability, as will be shown below. 

The collective dynamics throughout Stage II is very slow in the weak-coupling regime. 
Consider now (j6.1|) for two neighbouring clusters n and m which have both survived 
Stage I. The contribution of cluster m to the large parenthesis in the right-hand side 
of (|6.1|) is proportional to agx m . In the absence of coupling, we have x m ~ e^ 2 " -1 ** 5 / 2 , by 
virtue of (|3.4|) . The characteristic time scale of Stage II is reached when the product gx m 
becomes of order unity. It reads therefore 

? 1 

In-, (6.2) 



2a - 1 g 



t c ~t g-W«-V. (6.3) 

The separation of time scales between the fast individual and the slow collective dy- 
namics is therefore again parametrically large in the weak-coupling limit, although the 
divergence of the collective time scale is much less pronounced than in the mean-field 
limit [see f!5.19|) ]. The glassiness of the dynamics, with its manifest two-step relaxation, 
is illustrated in Figure El This figure shows a plot of the decay of the survival probability 
S(s) in one dimension. The dynamical equations (J6.1)) have been integrated numerically 
for a chain of 10 6 clusters, until every surviving cluster is isolated (see below). Both 
stages of the dynamics appear clearly on the plot, as well as the expected plateau value 
= 0.8, and the occurrence of a non-trivial limit survival probability Sr^ ~ 0.4134. 
Each curve corresponds to an interaction strength g a decade apart from its neighbour. 
It is accordingly shifted by 2 In 10 (thick bar), in accord with the estimate (J6.2|) . 

At the end of Stage II of the dynamics, i.e., in practice after a very long time, the 
system is left in a non-trivial attractor, which consists in a pattern where each cluster 
is isolated: all its first neighbours are dead, and is therefore a survivor, it survives and 
keeps growing forever. In the following, we shall call these attractors metastable states, 
in analogy with a variety of statistical-mechanical systems, where metastable states have 
been identified under various names in different contexts: valleys j^J, pure states [HI El, 
inherent structures [H], quasi-states 0. The common feature of metastable states in all 
these situations is that their number M generically grows exponentially with the system 
size (number of sites) N: 

W~exp(iV£). (6.4) 



17 



1.0 




2 1 1 1 1 1 1 1 1 1 1 1 

'0 10 20 30 40 50 

S 

Figure 3: Plot of the survival probability S(s) on the chain with Sm = 0.8. Left to right: 
Full line: g = 1CP 3 . Dashed line: g = 10 -4 . Long-dashed line: g = 10~ 5 . Dash-dotted line: 
g = 1CT 6 . The thick bar has length 2 In 10 = 4.605 (see text). 

The quantity £ is usually referred to as the configurational entropy, or complexity. 

The rest of this section is devoted to various characteristics of these attractors, such 
as their density (equal to the limit survival probability), spatial patterns, spatial corre- 
lations, and mass distribution of survivors. 

6.2 Limit survival probability 

The limit survival probability S'(oo), already emphasised in Figure El is just the density 
of a typical attractor, i.e., the fraction of the initial clusters which survive forever and 
take part in the attractor. The limit survival probability obeys the inequalities 

^(oo) < S(X), 5 (oo ) < 1/2. (6.5) 

The first inequality expresses that clusters can only disappear: the difference 1 — 5(i) 
(resp. 5(i) — 5(00)) is the fraction of clusters which die out during Stage I (resp. Stage II). 
The second inequality is a consequence of the fact that each surviving cluster is isolated. 
The densest configuration of lattice sites obeying this condition consists in occupying all 
the sites of either of the two sublattices, whose density is exactly 1/2. This value 1/2 of 
the highest density holds for the large family of so-called bipartite lattices, which includes 
hypercubic lattices (chain, square lattice, cubic lattice, ...). It is, however, not universal, 
and would e.g. be only 1/3 for the triangular lattice. 

In a given class of initial mass distributions, the limit survival probability <S7oo) is a 
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monotonically increasing function of the plateau value starting from S^) — for 
5(1) = 0, and going to a non-trivial maximum value S(oo)max < 1/2 in the Sn\ — > 1 limit. 

In the regime where Sn-\ is small, it can be shown that S^) is a l so small, and that it 
depends on alone. To do so, let us introduce the concept of supercluster. In analogy 
with a percolation cluster, a supercluster is defined as a set of k > 1 connected clusters 
which have survived Stage I, and such that all their neighbours have disappeared during 
Stage I. The fate of superclusters depends on their size k as follows. 

•k k = 1: If a supercluster consists of a single isolated cluster, it evolves in Stage II 
according to the dynamics of Section El it is a survivor, because its reduced square 
mass exceeds the threshold y* of ()3.5|) . For independent of initial masses, a super- 
cluster with k = 1 occurs with density p\ = S(i)(l — ) 2D . 

~k k — 2: If a supercluster consists of a pair of neighbouring clusters (represented 
as ••) both clusters evolve according to the dynamics of Section l4~2l the smaller dies 
out, while the larger is a survivor. We are thus left with «o or o« in the late stages 
of the dynamics. Such an event takes place with density p 2 = ~ S^) 2 ^ 20 ^. 

* k > 3: If three or more surviving clusters form a supercluster, they may a priori 
have more than one possible fate. Consider for instance a linear supercluster of 
three clusters (•••). If the middle one disappears first (•o«), the two end ones 
are isolated, and both will be survivors. If one of the end ones disappears first 
(e.g. ••o), the other two form an interacting pair, and only the larger of those two 
will survive forever (e.g. »oo). The pattern of the survivors, and even their number, 
therefore cannot be predicted a priori. 

The above enumeration implies S^) — p\ + p 2 /2 H , where the dots stand for the 

unknown contribution of superclusters with k > 3. As pi ~ S'(i), p 2 ~ S 2 ^, and so on, 
we are left with the expansion 

S {oD) = S (1) -DS 2 {1) + --- (6.6) 

The dependence of on details of the initial mass distribution at fixed therefore 
only appears at order Sf^. 

In the converse limit Sm — > 1, the limit survival probability reaches a non-trivial 
maximum value 5'( 00 ) max < 1/2, which depends very weakly on the mass distribution. 
For instance, in one dimension one has .S^ocOmax ~ 0.441 for an exponential distribution 
and S'( 0O ) max ~ 0.446 for a uniform distribution. Figure |U shows a plot of the limit survival 
probability SV^) (fraction of clusters that survive both fast and slow dynamics) against 
the plateau survival probability Sm (fraction of clusters that survive the fast dynamics 
of Stage I), for an exponential mass distribution in one, two, and three dimensions. 
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Figure 4: Plot of the limit survival probability against the plateau survival probabil- 

ity Sny Top to bottom: one dimension (circles), two dimensions (squares), three dimensions 
(triangles). Full lines (hardly visible through symbols): rational fits based on [2/2] Pade ap- 
proximants incorporating both terms of the expansion (|6.6|) . 



6.3 Spatial patterns of attractors 

We have seen already that the system generically ends up trapped in an attractor, or 
metastable state, where each cluster is isolated, i.e., surrounded by empty sites, and there- 
fore grows forever. The spatial patterns generated by survivors are therefore absolutely 
fixed once they are created by the dynamics. 

We now turn to a descriptive investigation of these patterns. One of our main ratio- 
nales for this exploration of spatial patterns derives from the cosmological model which 
was at the origin of the present model [Sj. In that context, it was of interest to obtain 
both the mass distribution of black holes and their spatial pattern. 

For ease of visualisation, we consider a square lattice. In the limit of highest density 
('S'(oo) — 1/2), there are only two possible 'ground-state' configurations of the system: one 
where the first sublattice is full of survivors, while the second is empty, and vice-versa. 
In this limit, patterns of surviving clusters are therefore perfect checkerboards. 

In order to describe the patterns of attractors, we are led to introduce at every site n, 
with integer co-ordinates (ni,ri2, . . . , no), both the survival index 



Or 



1 if the cluster at site n is a survivor, 
else, 



and the checkerboard index 



1) 



+niH \-n D 



(6.7) 



(6i 
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The survival index depicts very simply the pattern of surviving clusters surrounded by 
empty sites. The checkerboard index, on the other hand, represents, for each site, the 
local choice of one of the two symmetry-related 'ground states', i.e., of one of the two 
sublattices. This is easiest to understand using a one-dimensional example: the two 

ground states are H 1 1 or — I 1 ■■ All the 4> n are equal to —1 in the 

first pattern, and equal to +1 in the second pattern. The checkerboard index n thus 
classifies each site according to the particular ground state selected locally at this site. 

If the initial masses are large enough, so that the plateau survival probability Sm 
after Stage I is close to unity, the limit survival probability S(pa) is not far from its 
'ideal' highest value of 1/2. In this regime, attractors clearly exhibit a local checkerboard 
structure, as well as frozen-in defects with respect to a perfect checkerboard. The random 
structure of defects is entirely inherited from the random distribution of initial masses, 
because the dynamics is deterministic. 

Figure 03 shows a map of the survival index and of the checkerboard index for the 
same attractor of a 150 2 sample of the square lattice. This attractor has a density 
<3(oo) ~ 0.371. For greater clarity, we also zoom into a part of size 40 2 , in order to show 
better the correspondence between patterns of the survival and checkerboard indices. 
The frozen-in defects cause what appear to be little rivulets of voids which surround 
patches of perfect checkerboard. These islands of checkerboard are represented by black 
or white patches in the lower right-hand figure, depending on their parity, which is clearly 
visible from a comparison of the two figures. 

6.4 Spatial correlations 

The main consequence of short-range interactions is the generation of correlations be- 
tween clusters. In our model, a study of such correlations is especially meaningful in the 
long-time limit, when all the clusters which are still present are in fact survivors. In this 
regime, neighbouring sites are fully anticorrelated, because each survivor is surrounded 
by voids. However, at least close to the limit S'(oo) = 1/2, the next-nearest neighbours of 
a surviving cluster are expected to contain another survivor with high probability. Also, 
most survivors at late times should be quite massive, as they have both survived Stage I 
and then Stage II. We may thus expect that survival and mass correlations exhibit a 
rather similar dependence on the distance. These expectations are borne out by the 
following detailed study of correlation functions. 

Let us introduce the two-point correlation functions of the survival index and of the 
reduced mass at separation n: 



where (. . .) stands for a normalised ensemble average at a very late stage of the dynamics. 
These correlations are normalised so as to have C CT (0) = C x (0) = 1, whereas C a (e) = 




(6.9) 
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Figure 5: Two complementary representations of a typical pattern of surviving clusters on the 
square lattice, with Sn\ = 0.9, so that S^) ~ 0.371. Left: Map of the survival index. Black 
(resp. white) squares represent <r n = 1 (resp. cr n = 0), i.e., surviving (resp. dead) sites. Right: 
Map of the checkerboard index. Black (resp. white) squares represent n = +1 (resp. cp n = — 1). 
Top panels show a 150 2 sample. Bottom panels show an enlargement of a 40 2 region near the 
centre of the sample. 



C x (e) = 0, where e is any unit vector of the lattice. Their disconnected parts read 

^(oo) = lim C CT (n) = 5 (oo) , 0,(00) = , lim C x (n) = (6.10) 

n|^oo |n|^oo {X ) 

It is worth noticing that the correlations of the checkerboard index are not indepen- 
dent of those of the survival index. We have indeed: 

(0o0n) = (-l) ni+ - +nD (4^ (oo) a(n) - 4S (oo) + 1). (6.11) 

The existence of such an identity is quite natural, although it may seem surprising at 
first sight. Indeed the checkerboard index is only a different bookkeeping method for the 
same data on the positions of the survivors. 

Figure El shows a plot of the correlation functions C a (n) and C x (n) against distance, in 
one dimension. Both correlation functions exhibit a fast oscillatory convergence toward 
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Figure 6: Plot of the correlation functions against the distance n along the chain. Empty 
symbols: correlation C a {n) of the survival index. Full symbols: correlation C x {n) of the reduced 
mass. 

their disconnected parts. Similar features are observed for the on- axis correlations C a (ne) 
and C x {ne) on the square and cubic lattices. Figure [3 shows a plot of the logarithm of 
the absolute on-axis connected correlations C°(ne) = C a (ne) — CV(oo) and C^ine) = 
C x {ne) — C x (oo), against distance, in one, two, and three dimensions. The connected 
correlations are observed to fall off very fast to zero, so fast that it is hard to fit the 
precise form of their asymptotic decay. Neither a conventional exponential fall-off nor 
a more exotic super-exponential behaviour can be ruled out from the available data. A 
more accurate investigation of this point will form the subject of future investigations. It 
is worth recalling that, in the context of zero-temperature dynamics of Ising spin chains, 
the super-exponential fall-off of correlations in metastable states has been emphasised as 
a signature of the generation of a non-trivial measure on the space of attractors i- e -? 
loosely speaking, of the violation of Edwards' flatness hypothesis [TT] . 

6.5 Mass distribution of survivors 

To conclude this section, we examine the mass distribution of the survivors. In the late 
stages of the dynamics, when every cluster is isolated, its reduced mass grows according 
to Section 3, i.e., x ~ e^ 2a ~ 1 " >s ^ 2 . The mean reduced mass ((x)) of the surviving clusters 
therefore exactly follows the same growth law. Hence it is natural to measure cluster 
masses with respect to their mean, and to introduce the ratios 



x 



n 



m 



n 



(6.12) 
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Figure 7: Logarithmic plot of the absolute connected on-axis correlation functions, against 
distance n. Left: survival correlation. Right: mass correlation. Top to bottom: one dimension 
(circles), two dimensions (squares), three dimensions (triangles). 
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All the rescaled variables X n are independent of time in the very late stages of the dy- 
namics. They are therefore expected to have a well-defined limit probability distribution 
P(oo)(X). The first moment of this distribution is identically ((X)) = 1, whereas ()6.1U|) 
implies 

« (x 2 ) _ C g (oo) 

((*» 2 (00) (X)2 C x (00) 

We have measured the values ((X 2 )) « 1.73, ((X 2 )) « 1.93, and ((X 2 )) « 2.10, respectively 
in one, two, and three dimensions. Figure |H1 shows a plot of the whole rescaled distribution 
P(oo)(X) of the masses of survivors in these three cases. This distribution is observed to 
be both rather structureless and weakly dependent on dimensionality. 

The cosmological origins of the present model were our main motivation for looking at 
the mass distribution of survivors. In that context, survivors would be the descendants 
of primordial black holes, which would form part of the dark matter in the Universe. 
The above results suggest that the population of survivors is essentially given by a single 
mass scale, growing as the mass of a single isolated cluster, whereas the superimposed 
cluster-to-cluster fluctuations are described by the harmless distribution P^(X). This 
qualitative picture clearly ignores any cosmological details, but may still be of general 
interest from the viewpoint of using statistical-mechanical methods to probe such issues. 



7 Discussion 

In the above we have presented and investigated in detail the many facets of a novel 
and very rich model at the interface between non-linear dynamics and non-equilibrium 
statistical mechanics. Despite its origins in a rather exotic context, namely accretion 
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Figure 8: Plot of the limit distribution P^^X) of the rescaled mass X of surviving clusters. 
Circles: one dimension. Squares: two dimensions. Triangles: three dimensions. 



dynamics of black holes in a brane world j2HH]; the present model is of potential interest 
in many other situations, at least as far as qualitative features are concerned. 

As recalled in the beginning of the Introduction, the premise of equilibration generally 
holds in most physical instances. The rare exceptions to this principle include classical 
systems with long-range forces, with the noticeable example of the large-scale structure 
of the Universe. Similar instabilities, where tiny initial differences get amplified forever, 
are also met in other sciences, with one well-known example being the rich- get-richer 
principle in economics ^2]- It has been realised more recently that non-equilibrium 
statistical-mechanical models may exhibit a similar phenomenon, where a single micro- 
scopic state acquires a large population by virtue of a condensation transition, even in 
one dimension. The appearance of such a condensate can be viewed as a classical and 
non-equilibrium analogue of Bose-Einstein condensation The scenario of survival of 
the biggest arising from our model can be viewed as an extreme example of this instability, 
where the condensate ends up containing the entire mass. We reiterate that the model is 
non- conserving, in the sense that the final cluster becomes eventually more massive than 
all the earlier ones put together. The physical reason for this is that the interaction term 
derives from a radiation field in its cosmological incarnation which can be regarded 
as a mass reservoir. 

The present model also has many other specific features of interest in each of the 
geometries considered; we address each one in turn. 

First, for a finite assembly of coupled clusters, the model provides an interesting ex- 
ample of a deterministic dynamical system describing the evolution of competing agents. 
This class of problems has been studied at length in biophysics, one particularly well- 
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known example being the Lotka-Volterra system in population dynamics as described by 
predator-prey models ^3]. The present model exhibits a whole variety of types of trajec- 
tories in the transient regime, encoded in Phases I to IV of the phase diagram shown in 
Figure ^ Another distinguishing feature of our model is that only one survivor remains 
after a sufficiently long time, leading to the description of this model as a winner-takes- all 
type of model, despite the lack of a conservation law. 

Next, the thermodynamical mean-field limit of our model clearly shows many features 
of glassy behaviour. The evolution consists of two successive stages: a fast individual dy- 
namics in Stage I, followed by a slow collective dynamics in Stage II. In the weak-coupling 
regime, the characteristic reduced time scale s c of the slow dynamics diverges as 1/g 2 
[see (j5.19j) ]: the two time scales are well-separated in a way which is very reminiscent of 
the a and f3 relaxations observed in most glassy systems pEJ. A particularly interesting 
feature in this context is the universality of Stage II asymptotics, such as (J5.24)) or (|5.25|) . 
and the more unusual universality of the prefactor C, which only depends on the tail ex- 
ponent of the initial mass distribution. A further distinguishing characteristic is that the 
time separation between fast and slow dynamics in our model is simply given in terms 
of the coupling constant, and becomes parametrically large in the weak-coupling regime; 
this is true both within and beyond mean- field theory [see respectively 1)5.19)1 and 1)6.2)) ]. 
In most conventional glassy systems, the separation of time scales is governed by the 
appearance of a slowly growing length scale L(t), associated with some kind of ordering. 
By contrast, the glassiness in the present situation has dynamical origins: our model 
has features that are similar to driven systems, where time-scale separation arises either 
from a non-zero drift velocity V ^Hj or a non-zero shear rate 7 |16j . Such time-scale 
separations become parametrically large with the divergence (as V — > or 7 — > for 
those models, and as g — > in the present case) of the slow time scale. 

Beyond mean-field, e.g. on finite-dimensional lattices with nearest-neighbour cou- 
plings, the principle of survival of the biggest only applies locally; thus isolated clusters 
of dissimilar sizes are able to survive independently. At least qualitatively, this recalls 
local screening mechanisms in a variety of growth models, including cluster aggregation 
in suspensions [TZ]. In our model, however, the screening is extreme, in the sense that 
growing clusters, once isolated, are survivors, i.e., survive and keep on growing forever. A 
direct consequence of this is that the model exhibits both aging and metastability in a way 
that is qualitatively similar to what is observed in the mean-field limit, even though the 
time scale of the slow stage of the dynamics diverges less rapidly in the weak-coupling 
regime. Furthermore, the aging phase gets interrupted, as the system gets eventually 
trapped in a metastable state where a finite fraction of the entire lattice is occupied 
by isolated clusters which survive forever. If the clusters are initially large enough, the 
density is only slightly below 1/2, so that the spatial pattern of survivors has a 
local checkerboard structure. While similar checkerboard patterns have been observed in 
coupled map lattices [TBI EH] j our attractors are distinguished by their absolute stability: 
they are created irreversibly by a deterministic dynamics from the fluctuations in the 
initial distribution of their initial masses. Once created, they then survive forever. Many 
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questions regarding the statistics of the metastable states thus obtained remain open. It 
would be most interesting to know whether they are generated with a uniform measure a 
la Edwards ^T] in an appropriately defined ensemble, or with a highly non-trivial one, 
as suggested by recent investigations of one- dimensional spin models [TTH l2Hj. 
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